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Abstract 

The paper studies the quantum mechanical Coulomb problem on a 3-sphere. We present a special 
parametrization of the ellipto-spheroidal coordinate system suitable for the separation of variables. 
- - - , After quantization we get the explicit form of the spectrum and present an algebraic equation for 

■ the eigenvalues of the Runge-Lentz vector. We also present the wave functions expressed via Heun 

polynomials. 

(N 

^ ■ 1 Introduction 

' The Coulomb problem on the flat space is one of the two special systems of a particle moving in the 

conservative central field, whose symmetry group is larger then 0{N)\l\. This is a reason why, together 
with the oscillator, this system is of very high importance in mathematical physics. Both being textbook 
examples of the integrable systems they admit separation of variables in several coordinate systems. 

Due to the big number of hidden symmetries the Coulomb problem preserves its property of integra- 
bility after numerous deformations. Reducing the symmetry of the system, such deformations, however, 
can model some physically interesting and still integrable systems. Perhaps, the most known deformation 
of the Coulomb problem is an additional linear term in the potential, which represents a homogeneous 
electric field and in practice leads to the, so called. Stark effect. The Hamiltonian has the following form: 
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^\ ' In contrast to the more symmetric Coulomb Hamiltonian, which admits separation of variables in 

, spherical and parabolic coordinates, the Schrodinger (Hamilton- Jacobi) equation for this Hamiltonian 

' can be separated only in parabolic coordinates: 

y 

>— ^ ' S^— r + z, T] — r~z, tan ip = — (2) 
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Although this system is classically integrable, on the quantum level its exact solution is unknown and 
usually the last term in ^ is considered as a small perturbation. Assuming this, as usual, they use the 
solution of the non-perturbed system to calculate the deviation of the energy levels and the transition 
coefficients. It turns out (Stark effect), that all the first order transition coefiicients are equal to zero 
. and the deviations of the energy levels are proportional to the eigenvalue of the third component of the 
Runge-Lentz vector [2]. 

Because of the high importance of the Coulomb problem and its physical applications, the study of 
their various deformations is an important task. For example, a good question is: how does the space 
curvature affect the quantum mechanical systems? A very simple and important example of curved spaces 
is the sphere. In this paper we will discuss the system on the 3-dimensional sphere, which is fixed by the 
condition 

xl + xl+xl+xj^l, (3) 

with x^ being coordinates of the ambient 4-dimensional Euclidean space. 

In 1940 Schrodinger suggested an integrable spherical analogue of the Coulomb problem p]. He also 
found the spectrum of the system using his factorization method. Although this system was studied by 
many different authors (see [1] , [5] , [6] , [7] etc.) and shares many common properties with its flat ancestor, 
there are still many questions to be understood. For example, there were many attempts to solve the 
Schrodinger equation, but up to our knowledge there are no fully satisfactory results. 

It is known (see [8] and the references therein) that the Hamiltonian of the 3-dimensional spherical 
Coulomb system can be derived from a 4-dimensional oscillator performing a reduction procedure. In 
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[9] the authors used this technique to construct an integrable generahzation of the spherical Coulomb 
system with an additional term, which is the analogue of the term describing the homogeneous electric 
field in the flat case, i.e. a term which in the limit of the flat space reduces to the usual electric term. 
Typically, after such a reduction we obtain a system with the presence of a U{1) Dirac monopole (see 
[TD]). The potential of the system written in the coordinates x has the following form: 

U — ~-f^ + EeiXQXs, x^ = xl + x\ + x\ (4) 

The additional magnetic charge, however, does not affect the symmetry of the system (i.e. all the 
generators and their commutators are preserved) 

In the present work we will try to solve the Schrodinger equation for the Coulomb problem on the 
3-sphere in the generalized parabolic coordinates, which will allow us in the future to calculate the 
physically more interesting Stark effect. After the separation of variables in the Schrodinger equation for 
the wave function we obtain so called Heun equation (see, [I1],[I3])- It is the most general second order 
linear differential equation with four regular singularities. It usually appears in complicated nonlinear 
physical systems. For a review see [T4] . 

In comparison to the hypergeometric equation, which appears in the flat case, the Heun equation is 
much less studied. Namely, as we will see below, the orthogonality relations for the wave functions do 
not coincide with the known ones for the Heun polynomials. 

At the end of this section, let us mention that the described procedure almost without changes can 
be applied also for the case of a 5-sphere (see [T5]). 

The paper is organized as follows. 

In the Second section we describe the coordinate system suitable for our purpose. 
In the Third section we present the quantization of the discussed system. 

In the Last section we summarize the obtained results and discuss the difficulties, which do not allow 
us to calculate the deviation of the energy levels. 



2 The coordinate system 

For proceeding to the solution of the Schrodinger equation we need to choose a suitable coordinate 
system. By analogy with the flat case we look for some generalization of the parabolic coordinates in 
which the equation admits separation of variables. The form of the electric term in ([T]) written in parabolic 
coordinates ([2]) (ee;(C ~ "n)) the potential Q prompt us to write down the following relations: 



XqX^ = x:i\/ 1 - x^ ^ , x^l-x^^—;^, tan = — - , (1) 
where x is defined in 

It turns out that the coordinate system defined in such a way is orthogonal (the metrics in these 
coordinates is diagonal) 



I ^+'' .. n 

/ 4ai-e) " 

I (2) 

\ 



and is appropriate for our purposes. 

It is clear that both ^ and 77 run in the interval [—1,1] and, therefore, one can denote 



^ = sin(pi, ?7 = -sin.^2, (3) 

where (ip\,Lpi) are some new angular coordinates. We have purposely put the sign in the definition 
of (/32, in order to simplify things below. Avoiding it will only change the angle. 
Thus, we get 

. ¥'l+</'2 </?l-</?2 I : : I : : . , 

Xo = sm ^ , = cos ^ , xi = y — sm (/Si sm 952 cos 95, a;2 = y — smc^si sin(^2 sm(y5, (4j 

which automatically resolves (jSj and gives us the domain of (^1,2: 

¥'ie[0,7r], ^2e[-7r,0]. (5) 
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Remark In order to better understand the geometrical meaning of the angles </3i,2, let us for a moment 
fix the angle setting it, say, to 0. In this case these angles can define a chart on 5^ (this is possible due 
to the existence of the first Hopf map). (|H) gives a hint for a simple geometrical interpretation: ipi — if 2 
and If I + (p2 are the angles between the radius vector and the Xq and axis respectively. The first one 
is the azimuthal angle 9 of spherical coordinates. 

We call the coordinate system Q a "parabolic" one because of its vivid analogy with the flat one. 
However, a more correct term is the prolate elliptic coordinate system, which is in another parametrization 
discussed and used for the same reason in (T6j . 

The metrics in these coordinates has the following form: 

( sin yi— sin y2 n 
sin c^i 


At the end of this section, let us present the Hamiltonian of the non-perturbed MICZ-Coulomb system 
in the coordinates (^,77, (/?): 

Since the aim of this work is to present the solution of the Coulomb problem on the sphere, to 
simplify things we have set the magnetic charge to 0. The presence of the monopole, however can be 
easily recovered by a transparent change of parameters (see pTi ) . For completeness we also present the 
additional term, which is responsible for the presence of the magnetic charge [S]: 




sPv + s^ / 1 + 1 - Vl - 



(8) 



This is everything we need to know about the coordinate system, in order to proceed to the quanti- 
zation procedure. 

3 Schrodinger's equation 

In order to quantize the Hamiltonian ([7]) we need to replace its kinetic term with the corresponding 
Laplace operator 

—— (y^g''—) , =C, y^ = V, V^^V. 9'' 93k =51 5 = det(gy), (9) 



and consider the problem of the eigenvectors and eigenvalues for the obtained operator H: 

H%l) = Etp. (10) 

Just as the Hamilton- Jacobi equation, the Schrodinger's equation for this Hamiltonian obviously 
admits a separation of variables. Because of the axial symmetry, we look for the wave function ip in the 
form 

V'(C,'7,¥') = /(e)3('7)e""^ meN (11) 



and, after the separation of variables we find 



4?(i - e)f"{o + 4(1 - 2e)no + ^ + + /?) /(o = o 

47y(l - ^^)g"{r^) + 4(1 - 2r^^)g'{i^) + (-2^1^ + I|! + - /?) .9(7?) = 



(12) 



For proceeding further, we need to simplify these equations by performing a transformation of vari- 
ables. It was already mentioned that, although the variables ^ and 77 are more similar to the flat parabolic 
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coordinates, the set of variables {ipi, (/?2, ^) defined in ([3]) seems to be more natural for our problem. How- 
ever, as we will see below, it is even more suitable to enlarge our configuration space to the complex plane 
adding two additional conditions. Namely, instead of {ipi,ip2) we define two complex coordinates 

^i = e'^S Z2 = e^^^ ^ = -*4^- (13) 

2zi 2z2 

with two obvious conditions 

zizi = Z2Z2 = 1. (14) 
This substitution transforms the equations to the following form: 

z,{z! 1)/" + 2zff + (m^^ + E^+ + f ) / = 0, 

Z2{zl - l)g" + 2zlg' + (m^^ + e'-^ + + f ) 5 = 

The equations are completely identical and, therefore, for the time being we will drop the label 
indicating the variable z. Here we should notice that the transformation Z2 — > —22 in the equation for g 
changes only the sign of /3. Thus, after the quantization of the separation constant /3, it should satisfy 
the fact that, if (3^ is some quantized value, then so is —jS^. 

Let us look for the solution in the following form: 

/ = (;22 _ ^)^^i(l-V^+^)i7/ (16) 

Such a substitution leads us to an equation for HI 

HI" +{- + ^ + ^)hI'+( , "^^'^ HI . 0, (17) 
z-l z + lj \z{z - + z) J 



where 

r = 1 - VI + £^ + «7, 



e™, = \m\ + 1, a = if 



l + \m\ + ^l±E F^-VT+E+^ ^ b = I + - VT±ETE±VT±EEE (18) 



These quantities obey the additional condition a + fo+l^F + A-l-e. 

The equation ^Tll is called Heun's differential equation [T1],[T3]. Its solution is known as Heun's 
function and denoted as i?^ (— 1, g, a, 6, F, A, z). It is analytical in the disk \z\ < 1 and has Maclaurin 
expansion 



fc=0 



CkzK (19) 



The first parameter — 1 represents the position of the third singularity of the equation. In our case it is 
fixed to —1, and, therefore, in the future it will be omitted. 

Substituting the power series expansion in the equation (|17p for the coefficients we get: 

rCi + qC„ = 0, RkCk+i - qCk + PkCk-i = 0, (20) 

where 

Rk = -{k + l){k + T), Pk = {k-l + a)ik-l + b). (21) 

For the general case of the Heun equation the recursive formula is a bit more complicated. However, 
since in our case A = e and the third singularity is located in —1 (symmetric to 1), we get ([201) • 

The solution of this system can be written as follows. Let us define Xi = RiPi^i. Then, we have: 

Ck = E (-l)'^-(^o, X,,..., Xk-2)q''-^' (22) 

n Rr 

where £s{Xq, Xi, . . . , Xfc_2) is a polynomial of the power s of the variables (Xq, . . . , Xj,_2) such that it 
is a sum of all the monomials Xj-^Xj^ ■■■Xj^ (0 < ji < j2 < ■■■ < js 1^ k — 2) so that V p, g G 0, 1, s 
\jp — iq\ > l(if jp ~ jq or jp = + 1, thcu this monomial does not enter in the sum). 
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This solution (fT9|) is analytical only in the disk \z\ < 1. In order to make it analytical also on the 
boundary |z| = 1 and normalizable, we should cut down the series turning the Heun function to Heun's 
polynomials. In order to achieve this, we should impose 



Cn+i = C„+2 = 0, C„ ^ 0, n e N. 
For k = n + 1 from the (|^D|) we get 

Pn = 0; ^ a = —n. 
This immediately gives us the energy spectrum 



(23) 
(24) 



En = {n + \m\){n + \m\ + 2) - 



4(71- 



and 



r„,„, = -n - |?7i| - 



I'y 



2{\m\ + n+ 1) ' 
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2(\m\ 



1 



(25) 



(26) 



For k = n from (PD|) we get a condition for possible values of q. Namely, q should be a solution of 
with k = n + 1. An alternative, but equivalent requirement is that q should be an eigenvalue of the 
trix A = 5i+ijRi + 5j+i^iPj. 

The first requirement simplifies the coefficients (|18|) and we notice a symmetry 



Rk 



-{1 + k) 
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Thus, the matrix A*^"' takes the following form: 



/ 



V 







Ro 






Ri 



(|m| 



1) 



(27) 



-Ri 




0_ 

-Ro 



Rn-l 





(28) 



This matrix, obviously obey the following property: 



(29) 



where /- is the n x n anti-diagonal identity matrix. On the other hand, from the form of the matrix 
it is obvious that if X is an eigenvector of the matrix ^'^"^ then IX = —X. Taking into account 
that = l„xn, we can write the following sequence of identities: 

= qX = lA^'^hlX = qlX = -A^'^hX = qlX = A^"^^ X ^ -qX q = ~q, (30) 

and, therefore q is always imaginary and, if q is an eigenvector of A^"\ then so is —q. This means that 
(3 = ±2iq{see (HH)) is a good separation constant, in complete agreement with p^ . 

It is known (see.e.g.^) that both in the classical and quantum cases /3 represents the eigenvalue of 
the third component of the Runge-Lentz vector. Thus, as it was expected, solving the Schrodinger's 
equation by separating the variables in parabolic coordinates gives us the eigenvalues of three constants 
of motion: the Hamiltonian, the third component of the angular momentum, the third component of the 
Runge-Lentz vector. 

Hereinafter we will label q with indices n, m and k to indicate that it is the fc-th eigenvalue of the 
third component of the Runge-Lentz vector in the nth energy state. 
Finally, for the wave function we have 



(31) 
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where 

A^n — 1 eiyi — 1 ^nm 

fnkm = {Zl - {Zl + Zj^ ^ Hp{qnkrn;-n,bn,n,Tn7n,A,n;Zi), . , 

9nkm = {Z2 - 1) 2 {Z2 + 1) ^ ^ Hp{qnkm; -n, bnm, Tnm, ^m] Z2) , 

with Cnkm denoting the normahzation constant. 

Before we proceed further, let us notice several useful facts about the Heun's polynomials. Firstly, 
since according to A = e, it is obvious that the transformation z — —z in pT)) changes only the sign 
of g, and, therefore, one can state that if we choose Co(<z) = Co(— g), then 

Ck{q) = (— l)'^Cfe(— g), (with all the remaining constants coinciding). (33) 

and 

^P(.Qnkrn^ ^nrai ^nra-^ ^ni: ^) — Qnkm^ ^nm: ^nm: ^rm z). ('^4) 

On the other hand, for calculating the complex conjugate of the Heun's polynomial we notice that 



n n ^ n ^ 

i^^=^Cfc^fc = z-«5]a-fc^'■^7fz-"5]C^z^ C^ = ^a-fe. (35) 

fc=0 fe=0 fe=0 

If we take the complex conjugate of ((20|) and use ([ST)) we find that the coefficients satisfy exactly 
the same equation as Ck and Cq = Cg. Therefore: 

Co 

Cj, = -^Cn-k — Ck- (36) 
From this relation we can make two conclusions: 



Hp{qnkm, 'n, bnm J ^ nrm ^rm z) — Z Hp{c[nkrm n^bn^m^^ n-rriT t z) , CqCq — CnCn- (37) 

This is all we need to know, in order to proceed further. 



4 Discussion and results 

In this section we discuss the unsolved problems and summarize the results. 

So far, we have obtained the energy spectrum (j25p and an algebraic equation for the eigenvalues of 
the third component of the Runge-Lentz vector. The latter is determined from an algebraic equation 

Cn+l = (38) 

where C is defined in (|22p . This is a polynomial equation of power n. For a given principle quantum 
number n it has 2n for even n and 2n — 1 for odd n solutions. Thus, the number of its possible values 
coincides with that in the flat case. Considering the Stark effect on the flat space, we get that the 
first order deviations of the energy are proportional to the third component of the Runge-Lentz vector. 
Hence, it would be good to have an analytical expression for this quantity. Although the equation psp 
is a polynomial one, it has a very special form. Let us mention the similarity between (|22p and the 
expansion of a polynomial via elementary symmetric polynomials of its solutions 

n n 

Pn = l[{x-Z,)=Y,i-'^Ye,iZi,---,Z.^)x"-'. (39) 

i=l s=0 

Such an analogy can give us a hope that it is possible to explicitly express /3(or q) via the system 
parameters. 

We have also got the wave function (|3ip of the problem, up to a normalization constant. As always, 
it should be determined from the normalization integral: 

J i^nkmi^n'k'm'dV — 6nn'Skk'Smm' , (40) 
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where 



(41) 
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According to ([S]) the integration path of zi and Z2 is the upper and lower semicircles of zz = 1 respectively. 

Since the obtained wave functions are eigenfunctions of a Hermitian operator, it is clear that they are 
orthogonal, i.e. the condition (^0]) for \n ~ n'\ + \ k — k'\ + \m — m'\ > is satisfied automatically. 

There is, however, a point to be clarified. Due to the multiplier e'™''' in (I5T]) the orthogonality by m 
is obvious. Let us expand the l.h.s. of assuming that m — m': 

Ipnkmi'n'k'mdV ^ / p(zi , Z2 )fi'l (^1 ) W2 (zi ) Wi (^2) W2 (^2 )dzi ^^2 = 0, (42) 



where 

p(zi, Z2) = {zi ~ Z2) ((1 - zi)(l - z,)f--' ((1 + zi)(l + z,)y--' (ziZ2) '""'^^'"''"' -^ f 1 + — ) (43) 

\ Z\Zl ) 

and, for simplicity, we have defined 

Wi[z) ^ Hp{qnkm,-n,bnm,TnTn,^m,z) , W2{z) = iJp (q„'fc'm, -71, r„/™, A„, z) . (44) 

On the other hand, the known orthogonality relations are written for two polynomials with the same 
parameters (F, A, e, b) and different accessory parameters q and the power n. Namely, if we denote 

wi(z) = i/p(g„fc,-n, 6, r, A,z), W2(z) = ^/"p (g^'fe' , T, A, z) . (45) 

and 

Pstd{zi,Z2) = (Zi - Z2) ((1 - Zi)(l - Z2))'^"' ((1 + Zi)(l + Z2))'"' (^1^2)^"\ (46) 

then (|42)) is proved to be satisfied (see [12],[13]). Please, notice the absence of the last multiplier from 

gSl). 

Thus, it is clear that the orthogonality relations for the obtained wave functions do not coincide with 
the general ones. It is possible, since we have a very particular set of parameters F, A, e, h. 

It seems to be a matter of a different research to prove directly the orthogonality. This question is 
crucial also for the calculation of the Stark effect. Namely, the first order deviation of the energy levels 
and the transition coefficients are expressed via normalized wave functions as follows: 

/^Enkn'k' = j Ipn'k'vaAlj^nkmdV, (47) 

where A is the operator of the small perturbation. In our case it look as follows: 



1 z^-1 



2zi 2z2 

It was already mentioned, that in the flat only the coefficients with n = n' (energy deviations) are 
different from zero. It would be interesting to know whether this property survives in the case of spherical 
background space or here we have a qualitative difference. 

The mentioned problems do not seem to be unsolvable and will be studied in our future works. 
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